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ABSTRACT 

The metal-poor (Z ~ l/lOO-Z©) damped Lyman alpha system (DLA) at redshift 
z abs = 3.04984 in the z em ~ 3.030 QSO SDSS J1419+0829 has near-ideal properties 
for an accurate determination of the primordial abundance of deuterium, (D/H) p . We 
have analysed a high-quality spectrum of this object with software specifically designed 
to deduce the best fitting value of D/H and to assess comprehensively the random 
and systematic errors affecting this determination. We find (D/H)dla = (2.535 ± 
0.05) x 10~ 5 which in turn implies ^b,o^ 2 = 0.0223 ± 0.0009, in very good agreement 
with f2 b ,ofr 2 (CMB) = 0.0222 ± 0.0004 deduced from the angular power spectrum 
of the cosmic microwave background. If the value in this DLA is indeed the true 
(D/H) p produced by Big-Bang nucleosynthesis (BBN), there may be no need to invoke 
non-standard physics nor early astration of D to bring together fib. oh 2 (BBN) and 
fib.o^fCMB). The scatter between most of the reported values of (D/H) p in the 
literature may be due largely to unaccounted systematic errors and biases. Further 
progress in this area will require a homogeneous set of data comparable to those 
reported here and analysed in a self-consistent manner. Such an endeavour, while 
observationally demanding, has the potential of improving our understanding of BBN 
physics, including the relevant nuclear reactions, and the subsequent processing of 4 He 
and 7 Li through stars. 
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1 INTRODUCTION 

It is well known that the relative abundances of the light 
elements created in the first few minutes of cosmic history 
depend on the product of the present-day density of baryons 
in units of the critical density, fib.o, an d h 2 , where h is 
the present-day value of the Hubble parameter measured in 
units of 100 km s" 1 Mpc" 1 . When analysed in conjunction 
with the the power spectrum of temperature anisotropies 
of the cosmic background radiation (CMB), light element 
abundances can improve constraints on cosmological param- 
eters, primarily the spectral index of primordial fluctuations 
(Pettini et al. 2008), and on the effective number of light 
fermion species, A e ff (e.g. Simha & Steigman 2008; Nollett 
& Holder 2011). 



* Based on observations collected at the European Organisation 
for Astronomical Research in the Southern Hemisphere, Chile 
[VLT programme ID 085.A-0109(A)]. 
f email: pettini@ast.cam.ac.uk 



Among the elements created in Big Bang nucleosynthe- 
sis (BBN), three have been the focus of most observational 
efforts: 4 He, D and 7 Li. It is a lingering concern that D is 
the only one of the three whose primordial abundance is 
consistent, within the errors, with the value of Q\>fih 2 de- 
duced from the analysis of the CMB. The poor sensitivity 
of the primordial abundance of 4 He to fib.o^ 2 requires it to 
be measured with a precision which may be beyond what 
can be realistically achieved, at least until all the sources of 
systematic error are fully understood (Molaro 2008; Aver, 
Olive & Skillman 2012). On the other hand, a satisfactory 
explanation for the factor of ~ 3 discrepancy between the 
abundances of 7 Li in the oldest stars in our Galaxy and the 
value expected from BBN has yet to be found (Sbordone et 
al. 2010; Iocco 2012; Nissen & Schuster 2012 and references 
therein). 

Even D is not entirely without its problems. As pre- 
sciently pointed out by Adams (1976), its primordial abun- 
dance is most effectively deduced from the D/H ratio in 
metal-poor hydrogen clouds giving rise to absorption lines 
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Figure 1. Selected metal lines in the z abs = 3.04984 DLA in the QSO SDSS J1419+0829, reproduced from Cooke et al. (2011). In each 
panel, the black histogram is the observed spectrum and the red continuous line is the theoretical line profile fitted to the data. Vertical 
tick marks above the spectrum indicate the velocities of the two absorption components, with parameters listed in section [2] The j/-axis 
scale is residual intensity. The normalized quasar continuum and zero level are shown by the blue long-dashed and green dashed lines, 
respectively. 



in the spectra of high redshift quasars (QSOs). This is 
a difficult measurement, requiring not only the combined 
power of 8-10 m aperture optical telescopes and high res- 
olution echelle spectrographs, but also the right combina- 
tion of redshift, hydrogen column density, and internal ve- 
locity structure of the absorbing gas. Consequently, there 
are currently only about 10 QSO absorbers where the D I 
lines have been unambigously detected free from interfer- 
ence from unrelated absorption. The mean value so derived, 
(log(D/H) p ) = -4.556 ± 0.034 (Fumagalli et al. 2011) im- 
plies fi b ,o/i 2 (BBN) = 0.0213 ± 0.0012 which agrees within 
the errors with Sl bj0 /i 2 (CMB) = 0.0222 ± 0.0004 (Keisler et 
al. 2011). Less satisfactory is the fact that even this prime 
set of D/H measures exhibits a wider dispersion about the 
mean than expected on the basis of the quoted errors. While 
the most straightforward explanation of the scatter is that 
systematic errors affecting some, or all, of the D/H mea- 
surements have been underestimated, more radical ideas in- 
volving early destruction of D have also been proposed (e.g. 
Olive et al. 2012). 

Given the key role played by D in our understanding 
of BBN, it is important to identify and study more QSO 
absorbers with the optimum characteristics for the determi- 
nation of the primordial D abundance (D/H) p . This is the 
subject of the present paper. 



1.1 Metal-poor damped Lyman alpha systems 

Among the different classes of QSO absorbers, the so-called 
'damped Lyman alpha systems' (DLAs) are potentially the 
best targets for reliable measures of (D/H) p . DLAs are 
characterised by high column densities of neutral gas, with 
N(H i) > 2x 10 20 atoms cm" 2 (Wolfe et al. 2005), and appear 
to be associated with galaxies at early stages of chemical 
evolution (Pontzen et al. 2008; Peroux et al. 2012; Krogager 
et al. 2012). Typical metallicities at redshifts z — 2-4 are 
Zola ^ 1/10-1/100 Z Q (Ellison et al. 2012), but with a tail 
that extends to Zola < 1/1000 Z Q (Penprase et al. 2010; 
Cooke et al. 2011). 

The number of known DLAs with metallicities less 
than 1/100 of solar has increased in recent years thanks 
to large scale sky surveys (Noterdaeme et al. 2009), and it 
is these systems that have proved to be choice candidates 
for measuring (D/H) p for the following reasons: (i) the low 



metallicities imply negligible D astration (Romano et al. 
2006; Prodanovic et al. 2010), justifying the assumption that 
(D/H)dla = (D/H) ; (ii) the high Hi column densities re- 
sult in detectable optical depths of D I absorption in many 
lines in the Lyman series. This both increases the accuracy 
of the measurement of N(D i) compared with lower column 
density systems, where D I absorption is detected only in the 
strongest Lyman lines (e.g. Buries & Tytler 1998a, b), and 
reduces the likelihood of contamination by unrelated absorp- 
tion; (iii) the empirical correlation between metallicity and 
absorption line width (Murphy et al. 2007; Prochaska et al. 
2008) implies that the kinematics of Z £ 1/100 Z Q DLAs 
are normally quiescent, with the absorption lines extending 
over a narrow velocity interval (Cooke et al. 2011). Such 
kinematics greatly facilitate resolving the —82 km s _1 iso- 
tope shift of the D I Lyman lines from the neighbouring H I 
absorption. 

Thus, several of the very metal-poor DLAs in the recent 
survey by Cooke et al. (2011) are promising candidates for 
further observations targeting the high order D I and H I 
Lyman lines redshifted into the ground-base UV spectral 
region (A obs < 4000 A). In this paper we report the analysis 
of what we consider to be the most suitable case in the 
sample by Cooke et al. (2011) for determining the primordial 
abundance of D with a higher precision than is normally 
achievable. 

The paper is arranged as follows. In section [2]we briefly 
describe the observations (which have already been reported 
by Cooke et al. 2011) and summarise the most important 
properties of the DLA. Sections [3] and [4] are devoted to the 
derivation of the D/H ratio and its error, respectively. In 
section [5] we examine the cosmological implications which 
follow if the value of D /H deduced here is indeed the best 
measure of the primordial abundance of deuterium. We com- 
ment on earlier measurements of (D/H) p in section|6l before 
concluding in section [7] 



2 OBSERVATIONS AND DATA REDUCTION 

The spectrum of the g = 18.9 mag QSO SDSS J1419+0829 
at Zem — 3.030 shows a damped Lyman alpha system at 
Zabs = 3.04984; the small redshift difference between Lya 
emission and absorption results in an unusually 'clean' Lya 
absorption profile, with a red wing which is essentially free 
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Figure 2. Portion of the UVES spectrum of the QSO SDSS J1419+0829 (black), together with the model fit (red). The la error spectrum 
is shown in blue (near the zero level). Vertical dash lines mark the positions of QSO spectral features, as indicated. Light blue labels 
denote emission lines at z cm = 2.98576, green labels emission lines at z om = 3.04224, and red labels emission lines at z a t, s = 3.04954. 



from overlapping absorption by intergalactic gas (see Fig. 7 
of Cooke et al. 2011). Consequently, the column density of 
Hi is particularly well determined. 

Details of the acquisition and reduction of the spectrum 
of this QSO can be found in Cooke et al. (2011); here, we 
summarise the most important points. The spectrum was 
recorded with the Ultraviolet and Visual Echelle Spectro- 
graph (UVES; Dekker et al. 2000) mounted on UT 2 of the 
Very Large Telescope facility. The observations were carried 
out in service mode and consisted of a series of 2700 s ex- 
posures; the total integration time was 29 800 s. The echelle 
spectra were reduced with the standard UVES pipeline and 
then combined with the software package uves_poplefQ 
which merges individual echelle orders and maps the data 
onto a vacuum heliocentric wavelength scale. Flux calibra- 
tion was achieved by reference to the Sloan Digital Sky Sur- 
vey spectrum of the QSO. 

The final reduced spectrum covers the wavelength range 
A obs = 3710-6652 A (Ao = 916-1642 A in the rest frame of 
the DLA) with resolving power R ~ 40 000 (velocity reso- 
lution FWHM~ 7.5 km s" 1 ) sampled with ~ 3 pixels. The 
signal-to-noise ratio (S/N) of the data is relatively high con- 
sidering the magnitude of the source and the high spectral 
dispersion: S/N = 43 at 5000 A, and S/N ^ 20 over the 
entire DLA Lyman series (these are values per pixel in the 
continuum) . 

Cooke et al. (2011) deduced log JV(H l)/cm~ 2 = 20.40± 
0.03. Their analysis of 19 transitions of N I, O I, Si II, and Fe II 
showed that the metal lines consist of two absorption com- 
ponents of approximately equal optical depth, separated by 
Av = 13.8 km s _1 at z abs = 3.049649 and 3.049835 respec- 
tively. The Doppler parameters of the two components are 
b = 3.5 and 6.4 km s _1 respectively, where b = %/2cr, and a is 
the 1-D velocity dispersion of the absorbers projected along 

1 uvbs_popler is available from 

http://astronomy.swin.edu.au/~mmurphy/UVES_popler 



the line of sight. The metallicity of the DLA is approximately 
1/100 of solar: [O/H] = -1.92 ±0.05, [Si/H] = -2.08 ±0.03, 
and [Fe/H] = —2.33 ± 0.04 in the usual notation whereby 
[X/H] DLA = log (X/H) DLA - log (X/H) q; In Fig. □ we have 
reproduced examples of the metal lines in the DLA. 



3 MEASUREMENT OF D/H 

For the reasons outlined above, the spectral characteristics 
of this DLA are near-ideal for an accurate determination 
of the D/H ratio and a realistic assessment of its precision. 
To take advantage of this opportunity, we have performed a 
detailed spectral analysis in which we identify all of the pa- 
rameters which potentially affect the derivation of the D/H 
ratio from the UVES spectrum of J1419±0829, and then use 
a purpose-built code to solve for all of them simultaneously 
(via a \ 2 statistic) and thereby deduce the best fitting value 
of D/H in the DLA. Random and systematic errors affect- 
ing this determination were estimated using Monte Carlo 
techniques. We now discuss the whole procedure in detail. 

3.1 QSO continuum 

The absorption lines in the DLA are seen against the con- 
tinuum and emission line radiation from the QSO; the ex- 
act level of this flux as a function of wavelength affects di- 
rectly the measured line optical depths. Reconstruction of 
the QSO intrinsic spectrum is particularly important for the 
Hi Lya absorption line which: (i) is blended with the QSO 
Lya and Nv A1240 emission lines and (ii) affords the best 
constraint on iV(Hl). We approximated the QSO flux in 
the wavelength regions of interest by modelling the QSO in- 
trinsic spectrum with a power-law continuum on which are 
superimposed a number of emission lines (see Fig. [2]). The 
free parameters that were varied to arrive at the optimum 
fit (i.e. the fit giving the minimum value of \ 2 ) are: 0) t ne 
slope and (ii) scaling constant of the power-law continuum, 
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Figure 3. The Lya region in J1419+0829. Top panel: Observed 
spectrum in black and best-fitting model QSO spectrum in red. 
Middle panel: The normalized spectrum, obtained by dividing 
the observed spectrum by the model spectrum, is shown in black 
together with the best fitting damped Lya absorption profile 
(see section I3.2|l in red. The neutral hydrogen column density 
is log JV(Hl)/cm- 2 = 20.391 ± 0.008. Bottom panel: Expanded 
central portion of the middle panel shown in the rest frame of the 
^abs = 3.04984 DLA. In all three panels the lc error spectrum is 
shown in blue. 



and (iii) the redshifts, (iv) amplitudes and (v) widths of the 
emission lines. Multiplet wavelengths for the emission lines 
were taken from Morton (2003). 

We found that three sets of emission line redshifts were 
required to reproduce adequately the QSO emission line 
spectrum: a component at z cm = 3.04224 is seen in emission 
lines of ionised species (Cm, Civ, Niv, Nv, Sim, Tim) 
and in Lya; a lower redshift component at z om = 2.98576 
further contributes to high ionisation lines (Civ and Nv) 
and to Lya, while emission lines from low ionisation species 
(Cn, Ol, and Sin) are centred at z cm = 3.04954. The final 
fit is shown in red in Fig. [2] where the locations of the differ- 
ent emission line components are colour-coded. The region 
around SilvA1397 was not included in the fit because this 
emission line falls in a gap between the CCD chips on the 
UVES detector. Fig. [3] shows the region near the Lya emis- 
sion+absorption composite before and after division by the 
QSO flux. 

Further fine adjustments were applied to the continuum 
level in the proximity of the absorption lines to be analysed 
by fitting low order polynomials to spectral intervals deemed 
to be free of absorption. 



analysis of these features and ensuing abundance determi- 
nations have been reported by Cooke et al. (2011). Also 
covered by our spectrum is the entire Lyman series of the 
DLA, from Lya to the Lyman limit. Beyond Lyl4 (i.e. for 
transitions between the ground state and levels with princi- 
pal quantum number n > 15) the line wavelengths differ by 
AAo < 0.5 A and can no longer be resolved from one another. 
However, in 8 out of the 14 transitions from Lya to Lyl4, 
the D I absorption can be separated from the corresponding 
H I and is clear of other blends (see Fig. [4} . The /-values 
of the transitions in which D I is available for analysis range 
from 0.013940 (LyS) to 0.000469 (Lyl4); this factor of ~ 30 
spread in line strength provides a wide baseline on the curve 
of growth for an accurate measurement of N(Dl). 



3.2.1 Absorption components 

While only two absorption components are evident in the 
metal lines analysed by Cooke et al. (2011; see Fig. [l|, 
an additional, redshifted component was found to con- 
tribute to the Hi absorption. This component, separated 
by +12.7 km s" 1 from the redshift of the DLA, has only 
3.5 x 10 -4 the column density of the DLA, which explains 
why its presence cannot be discerned in the metal lines (re- 
call that in this DLA O/H ~ 5 x 10~ 6 and other metals 
are even less abundant). Even though its contribution to 
the D I absorption is minimal, this additional component is 
included in the D I lines by construction since, as explained 
below, our fitting procedure solves directly for the Dl/Hl 
ratio, rather than separately for A r (D i) and A r (H i) in indi- 
vidual components. 

The absorption from each component is characterised 
by four parameters: (i) redshift z, (ii) column density N, (iii) 
temperature T, and (iv) large-scale velocity dispersion. The 
last of these is assumed to be the same for all atoms and ions 
in the DLA, and is quantified by a turbulent Doppler param- 
eter bturb- On the other hand, a given temperature will corre- 
spond to different thermal Doppler parameters for elements 
of different masses according to the relation b\ = 2kT/m 
where k is the Boltzmann constant. The two contributions 
to the line widths add in quadrature: b\ ot = 6 2 urb + b%. The 
observation that the two components are unresolved in the 
D I lines while they are partially resolved in the metal lines 
(compare Figs. Q] and U) reflects the larger values of 6t for 
the lighter D compared to the metals, and leads to the con- 
clusion that in this DLA thermal broadening makes a non- 
negligible contribution to the line width. 

Thus, in fitting the absorption line spectrum of the 
DLA, we varied simultaneously z, N, T, and 6turb for each of 
the two main components of the H I, D I and metal lines. (For 
the third component seen only in Hi we varied z, N, and 
btot). The parameters defining the QSO continuum in the 
proximity of the Lya emission line were also varied simulta- 
neously to find the set of all of these values which minimises 
the x 2 statistic between observed and computed line profiles 
and the corresponding best-fitting value of D/H0 



3.2 Absorption lines 

Our UVES spectrum of J1419+0829 includes 19 metal ab- 
sorption lines associated with the z abs = 3.04984 DLA; the 



2 An additional free parameter which was allowed to vary in the 
X 2 minimisation procedure is the spectral resolution. The nominal 
value measured from the widths of emission lines from the Th- 
Ar lamp used for wavelength calibration is R = A/ A A ~ 40 000. 



© 2012 RAS, MNRAS 000, HE]] 



The primordial abundance of Deuterium 5 



1.0 
0.5 
0.0 

1.0 
0.5 
0.0 

1.0 

0.5 
0.0 



i 1 i 1 r~ 

Ji l _±y a i jj i_u 




i 1 i 1 i 1 i 1 r 

J lLy0 i jj uj 




1 1 1 1 1 1 1 1 1 
III III H I Ly9 


1.0 










0.5 


i , t- r*^ 


JK^-i — — -i 


0.0 



I 1 I 1 T" 

J ILy? i jj lu 





H li^* 5 I JJ LU 





-xiifW^ - i ; - i 1 



1.0 
0.5 
0.0 

1.0 
0.5 
0.0 



H I Lye 
i 


I 

JLL 


1 i 1 i 1 i 

LU 






H I Ly6 
i 


i 

JLL 


1 i 1 i 1 i 



i 1 i 1 i 1 i 1 i 

, , i_l JJ LLI _ HILyll 


i , r 


> 





1.0 

0.5 




1.0 


i 1 i 1 i 1 i 1 i 
HILy? I,, I,, 








0.5 




\ ft 




0.0 




— h 






-200 - 


100 


100 200 



i 1 i 1 i 1 i 1 i 
II | „H I Lyl4 




i 







1.0 
0.5 
0.0 



-200 -100 100 
Velocity Relative to z abs = 3.049840 (km s" 1 ) 



200 



Figure 4. Lyman series lines in the 2 a b a = 3.04984 DLA. Transitions to higher energy levels than Lyl4 are too closely spaced in 
wavelength to be resolved. The black histogram is the observed spectrum, while the red continuous line is the model fit to the absorption 
features (see section [3]2]| . Vertical tick marks above the spectrum indicate the three absorption components contributing to the Hi (red) 
and D I (green) absorption. For completeness, we also show (dashed red line) the model H I and D I absorption superposed on blended 
lines which were not used to constrain the model parameters. Dotted red lines are used for fits to features partly blended with Dl and 
Hi Lyman lines. In all panels, the y-axis scale is residual intensity. 



However, since the spectrograph slit may be illuminated slightly 
differently by the lamp and the QSO, we preferred to treat the 
resolution as a free parameter, rather than assuming that it is 
known perfectly. Reassuringly, the value which was found to fit 



It is worth stressing here that this procedure is some- 
what different from that adopted in some of the previous 



best the shapes of the absorption lines is R 
expected. 



39500 ± 700, as 
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a Fixed to be the same in all three components. 



b The third component is undetected in the metal and D I lines. 



derivations of the deuterium abundance. Specifically, in most 
previous studies of D/H in DLAs, the column densities of 
Hi and Dl were determined separately, the former being 
largely derived from the damping wings of the Lya line and 
the latter from the available D I lines in the Lyman series. 
The ratio of these two column densities then gives D/H. In 
the present case, by fitting simultaneously all of the Hi, Dl 
and metal lines, as well as the QSO flux from which the 
absorption takes place, we used all of the information avail- 
able in the DLA absorption spectrum to deduce directly the 
value of D/H (assumed to be the same in all three absorption 
components) which best fits the data. 

In order to fit the data as explained, we devel- 
oped custom-made software based on the same prin- 
ciples as vpfit, the software package most commonly 
used in the analysis of interstellar absorption lines (see 
://www. ast.cam.ac.uk/~rfc/vpfit. html). Our ABSORP- 
TION LINE fitting (alfit) code, written in the PYTHON 
programming language, uses a modified version of the mpfit 
package (Markwardt 2009)0 mpfit employs a Levenberg- 
Marquardt technique to derive the set of model parameters 
that minimises the x 2 statistic. We tested alfit extensively 
using fake data with known model parameters, and against 
vpfit in simpler cases than that described here, to ensure 
that its output is fully compatible with that of previous 
analyses which used vpfit. 

The parameters of the best fitting model are collected 
in Table [1] and the corresponding line profiles are shown 
superposed on the observed spectral features in Fig. [4] It is 
interesting to note that the two components that make up 
the DLA (components 1 and 2 in Table [T]) have tempera- 
tures T ~ 11000 and 10 000 K respectively. A comparably 
high temperature has been reported recently by Carswell 
et al. (2012) for a similarly metal-poor DLA, while in other 
cases somewhat lower temperatures have been inferred (Pet- 
tini et al. 2008; Cooke et al. 2012). Temperatures of a few 
thousand degrees for metal-poor DLAs are consistent with 
the lack of 21 cm (Srianand et al. 2012; Ellison et al. 2012) 
and molecular hydrogen (Petitjean et al. 2000) absorption 
in such systems. We are not able to derive a value of the 
temperature for component 3 because it is detected only in 
H I and our method relies on the comparison of the widths 
of absorption lines from elements of differing masses. 

3 MPFIT was originally written for the INTERACTIVE DATA LAN- 
GUAGE (IDL) environment and has recently been converted into 
PYTHON by Mark Rivers and Sergey Koposov. 
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Figure 5. Histograms showing the distribution of values of 
logDl/Hl in 2000 Monte Carlo random realisations of the spec- 
trum of J1419+0829, as described in section [4] 

The column densities listed in Table [T] differ by only 
~ 0.02-0.03 dex from the values reported by Cooke et al. 
(2011) whose simpler analysis did not differentiate between 
turbulent and thermal broadening of the absorption lines. 
For the total neutral hydrogen column density in the DLA 
we obtain log A/(H l)/cm -2 = 20.391 ± 0.008. The best- 
fitting value of D/H is log (D i/H i) = -4.601 ± 0.008. 



4 ERROR DETERMINATION 

An important aspect of the present analysis is an in-depth 
assessment of the errors, both random and systematic, which 
apply to the value of D/H deduced above. The first thing 
to point out is that the errors listed in Table Q] are purely 
random, reflecting the effect that the lcr error spectrum has 
on the derivation of the model parameters; these errors are 
just the diagonal terms from the covariance matrix. 

We tested the validity of this estimate of the random 
error applicable to the best fitting value of D /H using a suite 
of simulations, as follows. The parameters defining the best 
fitting model were used to generate a fake UVES spectrum 
which was then perturbed by a random realisation of the 
error spectrum. To ensure that our choice of starting param- 
eters in the line fitting procedure does not bias the result, 
for each realisation we randomly drew a new set of start- 
ing parameters from the covariance array of the best-fitting 
model that describes the real data. The resulting spectrum 
was then analysed in an identical manner to the real spec- 



http 



© 2012 RAS, MNRAS 000,[2fTT] 



The primordial abundance of Deuterium 7 



Component 1 



„ -4.58 
SE 

a -4.60 



00 
O 



-4.62 



1 , 1 , 1 1 
Jog N(H 1) = 20.231_ 

; (•) ; 


z ab8 = 3.049840 

; (J) ; 


i ' i ' i 

- biui* = km/s _ 


_T D[A = l.l_3e + 04 K_ 

" ^9) ■ 


i i ~i , i 


i.i.i 


1,1,1 


1 I [ I 1 



-4.58 



- -4.60 



-4.62 



-0.02 0.02 -4xl0" 8 4xl0" 6 -0.4 0.4 -600 600 
Component 2 



-4.58 



-4.60 



GO 

o 



-4.62 - 



-, 1 1 1 , 1 r 

.log N(H I) = 19.86 





n — 1 — i — 1 — i — 

b turb = 23 km / s 





J I I 



-4.58 



-4.60 



- -4.62 



-0.04 0.04 -4xl0- 6 4xl0" 6 -0.4 0.4 -300 300 

Figure 6. Elliptical contour plots illustrating the covariance between model parameters and the Dl/Hl ratio. The black cross near the 
centre of each panel shows the best- fitting value of each parameter, as indicated (see also Table [TJ, while the ir-axis gives the deviation 
from this value. The red-shaded contours refer to the 1 and 2(7 ranges for the random errors; the corresponding ranges for the systematic 
errors are shown in blue. 



trum (as described in section [3]), to derive a new set of best 
fitting parameters, including the D/H value. This procedure 
was repeated 2000 times and the distribution of resulting 
values examined to infer the la random error on D/H (i.e. 
the range about the mean encompassing 68% of the 2000 
returned values). We found that the error on log(Dl/Hl) 
so derived is lcr ran = ±0.008 (see Fig. [5]), in good agree- 
ment with the value obtained from the diagonal terms of 
the covariance matrix. 

The two major sources of systematic error are the con- 
tinuum placement and the zero level; both affect the mea- 
sured optical depths of the absorption lines. To assess their 
impact on the derived value of D/H, we used the above 2000 
Monte Carlo realisations and applied continuum and zero 
level adjustments that reflect the operations performed on 
the real data. In this way, we can also identify any systematic 
bias that may be introduced as a result of our continuum 
and zero level placements. 

Examination of the black cores of saturated H I lines of 
the DLA's Lyman series indicated that the zero level: (i) 
is independent of wavelength, and (ii) has an uncertainty 
of ±0.0013 (Iff) of the continuum flux. For the 2000 real- 
isations described above, we perturbed the zero level by a 
constant amount, drawn at random from a Gaussian distri- 
bution of values about zero with a dispersion of 0.0013. The 
uncertainty resulting from the continuum placement was es- 
timated by repeating the same fine continuum adjustments 
described in section 13.11 on each of the 2000 randomly gen- 
erated spectra and redetermining the best fitting model pa- 
rameters as before. 

For each realisation, we take the difference between the 
best-fitting value of D/H, respectively with and without the 
aforementioned corrections, as an indication of the system- 



atic uncertainty affecting the derivation of the D/H ratio. 
The distribution of such differences in the 2000 trials then 
gives the estimate of the systematic error and bias in our 
measure of log (D i/H i) . 

As can be seen from Fig. O this second set of tests in- 
dicates that: (i) the magnitude of the systematic error is 
small: l<T sys = ±0.003; and (ii) there appears to be a bias 
of —0.005 in our estimate of log(Dl/Hl) introduced by a 
(small) error in the continuum placement. A further advan- 
tage of our approach is that it allows us to examine in detail 
the sensitivity of the D/H ratio to each of the model param- 
eters entering in its derivation. In Fig. [5] we have reproduced 
some examples; the error contours are 'well-behaved' and we 
see only mild correlations between Dl/Hl and some of the 
model parameters. 

The result of all of the above tests is as follows: 

log (D i/H i) = -4.596 ± 0.008 ± 0.003 (1) 

(random and systematic lcr errors respectively), which in- 
cludes the correction of ±0.005 dex for the bias in the con- 
tinuum and zero level placement. 

Combining the two errors in quadrature, leads to: 

log (D i/H i) = -4.596 ±0.009, (2) 

or: 

TO 5 (Di/Hi) = 2.535 ±0.05. (3) 



5 COSMOLOGICAL IMPLICATIONS 

Our determination of D/H in J1419±0829 is more precise 
than those reported up to now in other QSO absorbers for 
the following reasons: (i) The red wing of the damped Lya 
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line suffers very little contamination by gas unrelated to the 
DLA, easing the determination of iV(Hl); (ii) eight Dl Ly- 
man lines of widely differing /-values are accessible; (iii) the 
kinematic structure of the gas is simple, with only 2-3 com- 
ponents contributing to the absorption lines; and (iv) the 
spectrum analysed is of moderately high S/N. We therefore 
consider it worthwhile to examine the cosmological impli- 
cations of the new measurement reported here before dis- 
cussing the full sample of available D/H measures at high 
redshift. 

In the following, we take: 



(Di/Hi) 



(D/H) DLA = (D/H) p . 



(4) 



The assumptions underlying these equalities are that: (i) the 
fractional ionizations of H and D are the same, (ii) D is not 
depleted relative to H, and (iii) the destruction of D through 
astration prior to the time when we observe the DLA has 
been negligible. Concerning the first assumption, we are not 
aware of a physical process that would under- or over-ionise 
one isotope relative to the other. Dust depletion of D in 
the local interstellar medium has been proposed to explain 
the surprising range of D/H values found along different 
sightlines in our Galaxy (Linsky et al. 2006) , but is unlikely 
to be important in metal- and dust-poor DLAs where even 
highly refractory elements are present in near-solar relative 
proportions (Akerman et al. 2005; Vladilo et al. 2006; Ellison 
et al. 2007). Given the low metallicity of the z abs = 3.04984 
DLA, where N, O, Si, and Fe have abundances less than 
~ 1/100 of solar (Cooke et al. 2011), the third assumption 
is supported by chemical evolution models which entertain 
little reduction of the D abundance from its primordial value 
when such a small fraction of the gas has evidently been 
cycled through stars (Romano et al. 2006). 

Recently, Steigman (private communication) has up- 
dated the relations between (D/H) p and Qb,oh 2 given by 
Simha & Steigman (2008) and Steigman (2007), as follows: 



10 5 (D/H) = 2.60(1 ±0.06) 



ma - 6(5 - 1) 



(5) 



where tjio is the post-e annihilation ratio of the numbers 
of baryons and photons in units of 10 -10 : 

r/io = 10 10 n b /n 7 = 273.9 fi b ,o h 2 , (6) 

and 5 is the non-standard expansion rate factor, related 
to the number of additional (equivalent) neutrinos AN„ = 
N„ - 3 by: 

2 



1 + 



7AN„ 
43 



The 6% error in the conversion factor in eq. ([5]) reflects the 
uncertainties in the nuclear reaction rates, particularly the 
d(p,7) 3 He cross-section (e.g. Nollett & Holder 2011), used 
in BBN codes (Steigman, private communication). 

For standard BBN with 5=1, eqs. ©, J5j, and @ 
yield: 



1000 b ,o/i (BBN) = 2.23 ± 0.03 ± 0.08 



(8) 



where the error terms reflect the uncertainties in, respec- 
tively, (D/H) p (eq. [3]) and BBN calculations (eq. [5]). Com- 
bining the two error terms in quadrature, we have: 
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Figure 7. Measures of the deuterium abundance in high red- 
shift QSO absorbers. Only cases were the deuterium absorption 
is clearly resolved from nearby spectral features are shown (see 
text). The red star refers to the new measurement reported here, 
with errors smaller than the symbol size. The horizontal lines are 
drawn at the weighted mean value of log (D/H) and its error, as 
determined with the bootstrap method. The yellow shaded area 
shows the range in Q b /i 2 (CMB) from Keisler et al. (2011). 



Recently, Keisler et al. (2011) combined their measure- 
ment of the CMB angular power spectrum from the South 
Pole Telescope (SPT) with the power spectra from the seven- 
year Wilkinson Microwave Anisotropy Probe (WMAP) data 
release to better constrain cosmological parameters. From 
this analysis, it was concluded: 



100fi b ,o/i (CMB) = 2.22 ± 0.042 



(10) 



(see Table 3 of Keisler et al. 2011). 

The agreement between eqs. ([9} and (|10p is very en- 
couraging. If the value of (D/H) p we have deduced here is 
indeed the correct one, there may be no need to appeal 
to non-standard physics to reconcile BBN and CMB esti- 
mates of flb,o (see, for example, the discussion of this point 
by Hamann et al. 2011). Specifically, if we adopt the value 
100n b ,o/i 2 (CMB) = 2.22 ± 0.042 from Keisler et al. (2011) 
and then use eqs. ©, ©, and to solve for AN„ using 
the determination of (D/H) p reported here (eq. [3jl , we find: 



N„ = 3.0 ±0.5 



(11) 



lOOfib.oft (BBN) = 2.23 ± 0.09 . 



(9) 



combining in quadrature the errors on (D/H) p , BBN reac- 
tion rates and f2 b; o/i 2 (CMB), with the BBN rates making 
the largest contribution to the total uncertainty. This value 
is in good agreement with N„ = 2.984 ±0.008 deduced from 
the width of Z-boson decays in electron-positron colliders 
(ALEPH Collaboration et al. 2006). 



6 COMPARISON WITH OTHER 
MEASUREMENTS OF D/H 

In the previous section, we used our new, precise measure 
of D/H to compare the present day baryon density derived 
from BBN calculations and CMB measurements. However, a 
full comparison between flb,o(BBN) and S7 bj o(CMB) should 
include all reliable measurements of the primordial abun- 
dances of the light elements created in BBN. It is beyond 
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Table 2. Prime Sample of D/H measurements in QSO Absorption Line Systems 









lOti i V 1 11 1 t 

( pm ^\ 
l Lm ) 




log yu/ nj 


Ref. b 


HS 0105+1619 


2.640 


2.53600 


19.42 ± 0.01 


-1.73 


—4.60 ± 0.04 


1 


Q0913+072 


2.785 


2.61843 


20.34 ± 0.04 


-2.40 


-4.56 ±0.04 


2, 3 


Q1009+299 


2.640 


2.50357 


17.39 ±0.06 


< -0.70 c 


-4.40 ±0.07 


4 


SDSS J1134+5742 


3.522 


3.41088 


17.95 ±0.05 


< -1.9 d 


-4.69 ±0.13 


5 


Q1243+307 


2.558 


2.52566 


19.73 ±0.04 


-2.79 


-4.62 ±0.05 


6 


SDSS J1337+3152 


3.174 


3.16768 


20.41 ±0.15 


-2.68 


-4.93 ±0.15 


7 


SDSS J1419+0829 


3.030 


3.04984 


20.391 ±0.008 


-1.92 


-4.596 ± 0.009 


8, 9 


SDSS J1558-0031 


2.823 


2.70262 


20.67 ± 0.05 


-1.50 


-4.48 ± 0.06 


10 


Q1937-101 


3.787 


3.25601 


18.26 ±0.02 


-2.0° 


-4.80 ±0.08 


11 


Q1937-101 


3.787 


3.57220 


17.86 ±0.02 


< -0.9 


-4.48 ± 0.04 


12 


Q2206-199 


2.559 


2.07624 


20.43 ±0.04 


-2.07 


-4.78 ±0.09 


2, 13 



a Rclative to the solar value log(O/H) + 12 = 8.69 (Asplund et al. 2009). 

b References - (1) O'Meara et al. (2001), (2) Pettini et al. (2008a), (3) Pettini et al. (2008b), 
(4) Buries & Tytler (1998b), (5) Fumagalli et al. (2011), (6) Kirkman et al. (2003), (7) Srianand 
et al. (2010), (8) This work, (9) Cooke et al. (2011), (10) O'Meara et al. (2006), (11) Crighton 
et al. (2004), (12) Buries & Tytler (1998a), (13) Pettini & Bowen (2001). 

c This is a very conservative upper limit on the metallicity. Buries & Tytler (1998b) estimate 
[Si/H] ~ —2.5 and [C/H] ~ —2.9 from photoionisation modelling. 

d This is a conservative upper limit on the metallicity. Fumagalli et al. (2011) estimate [Si/H] ~ 
—4.2 from photoionisation modelling. 

°This value refers to [Si/H] as the O abundance was not measured by Crighton et al. (2004). 



the scope of this paper to comment on the difficulties in 
deducing the primordial abundances of 4 He and 7 Li from, 
respectively, the emission line spectra of metal-poor H n re- 
gions and the absorption spectra of some of the oldest stars 
in the Galaxy. We refer the interested reader to the review 
by Steigman (2007) for a detailed discussion. 

Here, we focus specifically on the derivation of (D/H) p . 
Since the compilation assembled by Pettini et al. (2008), 
two new measurements of (D/H) p have been reported, by 
Srianand et al. (2010) and Fumagalli et al. (2011). In the 
former case, the Dl lines are only partially resolved from 
nearby blends, so that its inclusion in what we consider to 
be the most reliable set of measurements is questionable but, 
for completeness, we have included this DLA in the sample. 
Together with the new detection reported here, the full data 
set now consists of 11 measurements of (D/H) p whose main 
characteristics are collected in Table [2] 

As can be seen from Fig. [7] there is a troublesome dis- 
persion between the 11 measurements, well in excess of the 
quoted errors. The problem is that even this set is highly het- 
erogeneous. The number of Dl Lyman lines covered varies 
between the 11 QSO absorbers, with only Lya and Ly,9 
available in systems of lower A r (Hl). The degree of line 
blending is significantly worse in some cases than in oth- 
ers. Some spectra are of inferior S/N ratio than others, such 
as that of the z abs = 2.07624 DLA in Q2206-199 which, be- 
cause of the lower redshift of this absorber, necessitated ob- 
servations at ultraviolet wavelengths with the Hubble Space 
Telescope (Pettini & Bowen 2001). Perhaps most impor- 
tantly, the efforts devoted to evaluating realistic errors vary 
considerably between the different analyses collected in Ta- 
ble O and we suspect that in at least some cases the errors 
quoted may well be underestimates, and the values reported 
may suffer from biases that are unaccounted for. 



It is our contention that the only way forward with 
this problem is to assemble a uniform sample of metal-poor 
DLAs, observed at comparable S/N, and analysed consis- 
tently following procedures similar to those employed here. 
Unfortunately, this ambitious goal is still some way off in the 
future. Given the reservations expressed above, it is unclear 
to the present authors whether there is much to be learnt 
by averaging the measurements in Table [2] Nevertheless, for 
completeness, we have estimated the mean ((D/H) p ) using 
the bootstrap method. That is, we calculated the weighted 
mean of the 11 measurements in Table [2] adopting the errors 
quoted in the original reports, and repeated the calculation 
10000 times by random sampling (with substitution) of the 
11 values of (D/H) p . The mean and standard deviation of 
the results of 10000 such trials are: 

(log (D/H) p ) = -4.585 ±0.02. (12) 

For comparison, Fumagalli et al. (2011) calculated 
(log(D/H) p ) = -4.556 ± 0.034, while Pettini et al. (2008) 
reported (log(D/H) p ) = —4.55 ± 0.03; both estimates are 
from subsets of the data in Table [2] Thus, the new case 
highlighted here has reduced slightly both the mean of the 
sample and its error, but the main conclusion is that the 
value of (log(D/H) p ) has remained relatively stable over the 
last few years. 



7 SUMMARY AND CONCLUSIONS 

In summary, we have identified a metal-poor DLA from 
the recent survey of such systems by Cooke et al. (2011) 
with near-ideal properties for an accurate determination of 
the primordial abundance of deuterium. To capitalise on 
this rare opportunity, we have developed a spectral anal- 
ysis specifically targeted at finding the most likely value of 
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D/H by varying all of the spectral parameters which have a 
bearing on this measurement. We have also paid particular 
attention to the assessment of the errors and biases, both 
random and systematic, which affect our determination of 
D/H. 

Our principal result is that the value of (D/H) p so de- 
rived, (D/H) p = (2.535 ±0.05) x 10~ 5 , implies Q, hfi h 2 = 
0.0223 ± 0.0009 which is in excellent agreement with 
fib,oft 2 (CMB) = 0.0222 ± 0.0004 deduced from the joint 
analysis of the CMB angular power spectrum recorded with 
the WMAP and SPT experiments. On the basis of the new 
measurement reported here, there is no need to invoke non- 
standard physics, nor appeal to early astration of deuterium, 
to reconcile CMB and D/H determinations of the cosmic 
density of baryons. 

When considered in the context of previous reports of 
D/H values in QSO absorbers, the new case presented in this 
paper highlights the heterogeneity of the existing data set. 
The weighted mean of what are often considered to be the 
best 11 available measurements is ((D/H) ) = (2.6 ±0.1) x 
10~ 5 , but it is hard to assess the significance of this mean 
value when it is likely that systematic errors and biases may 
have been overlooked in at least some earlier analyses. 

Looking ahead, we can expect a further refinement in 
the determination of fib, o^ 2 (CMB) from the Planck mission 
(Planck Collaboration et al. 2011). However, improvements 
in the value of (D/H) p made possible by the identification 
and careful analysis of suitable metal-poor DLAs remain a 
priority, given their potential for clarifying our understand- 
ing of primordial nucleosynthesis (including the relevant nu- 
clear reaction rates) and of the subsequent processing of 4 He 
and 7 Li. As an added bonus, more detections of D I absorp- 
tion in metal-poor DLAs will increase the limited statistics 
on the kinetic temperature of these structures and help us 
isolate the dominant physical processes that regulate it. 



ACKNOWLEDGEMENTS 

We are grateful to the ESO time assignment committee for 
their continuing support of this demanding observational 
programme, and to the staff astronomers at Paranal who 
conducted the observations. It is a pleasure to acknowledge 
helpful and stimulating conversations with Bob Carswell, 
Dawn Erb, Jan Hamann, Paul Hewett, Mike Irwin, Lloyd 
Knox, Michael Murphy, Subir Sarkar, and Gary Steigman. 
We are grateful to Gary Steigman for communicating his re- 
vised fitting formulae in advance of publication, and to the 
referee, John Webb, for valuable comments which improved 
the paper. 



REFERENCES 

Adams T. F. 1976, A&A, 50, 461 

Akerman, C. J. Ellison S. L., Pettini M., Steidel C. C. 2005, 
A&A, 440, 499 

ALEPH Collaboration, DELPHI Collaboration, L3 Collabora- 
tion, et al. 2006, Phys. Rep., 427, 257 

Asplund M., Grevesse N., Sauval A. J., Scott, P., 2009, 
ARA&A, 47, 481 

Aver E., Olive K. A., Skillman E. D. 2012, JCAP, 4, 4 

Buries S., Tytler D. 1998a, ApJ, 499, 699 



Buries S., Tytler D. 1998b, ApJ, 507, 732 

Carswell R. F., Becker G. D., Jorgenson R. A., Murphy M. T., 

Wolfe A. M. 2012, MNRAS, 422, 1700 
Cooke R., Pettini M., Steidel C. C, Rudie G. C, Nissen P. E. 

2011, MNRAS, 417, 1534 

Cooke R., Pettini M., Murphy M. 2012, MNRAS, in press 

larXiv:1201.1004t 
Dckker H., D'Odorico S., Kaufer A., Delabre B., Kotzlowski 

H., 2000, SPIE, 4008, 534 
Ellison S. L., Prochaska J. X., Lopez S. 2007, MNRAS, 380, 

1245 

Ellison S. L., Kanekar N., Prochaska J. X., Momjian E., 

Worseck G., 2012, MNRAS, submitted 
Fumagalli M., O'Meara J. M., Prochaska, J.X. 2011, Science, 

334, 1245 

Hamann J., Hannestad S., Raffelt, G. G. Wong Y. Y. Y. 2011, 
JCAP, 9, 34 

Iocco, F. 2012, Mem. S.A.It., in press JarXiv:1206.2396ll 
Keisler R., Reichardt C. L., Aird K. A., et al. 2011, ApJ, 743, 
28 

Kirkman D., Tytler D, Suzuki N., O'Meara J. M., Lubin D. 

2003, ApJS, 149, 1 
Krogager J.-K., Fynbo J. P. U., M0ller P., et al. 2012, MNRAS, 

in press l larXiv:1204.2833H 
Linsky J. L., Draine B. T., Moos H. W., et al. 2006, ApJ, 647, 

1106 

Markwardt C. B., 2009, ASPC, 411, 251 

Molaro, P. 2008, in Knapen J. H., Mahoney T. J., Vazdekis A., 
eds, ASP Conf. Ser. Vol. 390, Pathways Through an Eclectic 
Universe, Astron. Soc. Pac, San Francisco, p. 472 

Morton, D. C. 2003, ApJS, 149, 205 

Murphy M. T., Curran S. J., Webb J. K., Menager H., Zych 

B. J. 2007, MNRAS, 376, 673 
Nissen P. E., Schuster W. J. 2012, A&A, 543, A28 
Nollett K. M., Holder G. P. 2011, Phys. Rev. D, submitted 

l|arXiv:1112.2683|l 
Noterdaeme P., Petitjean P., Ledoux C, Srianand R. 2009, 

A&A, 505, 1087 
Olive K. A., Petitjean P., Vangioni E., Silk J. 2012, 

larXiv:1203.5701l 

O'Meara J. M., Buries S., Prochaska J.X., Prochter G. E. 2006, 
ApJ, 649, L61 

O'Meara J. M., Tytler D., Kirkman D., Suzuki N., Prochaska 
J.X., Lubin D., Wolfe A. M. 2001, ApJ, 552, 718 

Penprase B. E., Prochaska J. X., Sargent W. L. W., Toro- 
Martinez I., Beeler D. J. 2010, ApJ, 721, 1 

Peroux C, Bouche N., Kulkarni V. P., York D. G., Vladilo G. 

2012, MNRAS, 419, 3060 

Petitjean P., Srianand R., & Ledoux C. 2000, A&A, 364, L26 

Pettini M., Bowen D. V. 2001, ApJ, 560, 41 

Pettini M., Zych B. J., Steidel C. C, Chaffee F. H. 2008a, 

MNRAS, 385, 2011 
Pettini M., Zych B. J, Murphy M. T., Lewis A., Steidel C. C. 

2008b, MNRAS, 391, 1499 
Planck Collaboration, Ade P. A. R., Aghanim N., et al. 2011, 

A&A, 536, Al 

Pontzen A., Governato F., Pettini M., et al. 2008, MNRAS, 
390, 1349 

Prochaska J. X., Chen H.-W., Wolfe A. M., Dessauges- 

Zavadsky M., Bloom J. S. 2008, ApJ, 672, 59 
Prodanovic T., Steigman G, Fields B. D. 2010, MNRAS, 406, 

1108 

Romano D., Tosi M., Chiappini C, Matteucci F. 2006, MN- 
RAS, 369, 295 

Sbordone L., Bonifacio P., Caffau E., et al. 2010, A&A, 522, 
A26 

Simha V., & Steigman G. 2008, JCAP, 6, 16 

Srianand R., Gupta N., Petitjean P., Noterdaeme P., Ledoux, 



© 2012 RAS, MNRAS 000,[THIT] 



The primordial abundance of Deuterium 

C. 2010, MNRAS, 405, 1888 
Srianand R., Gupta N., Petitjean P., et al. 2012, MNRAS, 421, 
651 

Steigman G. 2007, Annual Review of Nuclear and Particle Sci- 
ence, 57, 463 

Vladilo G., Centurion M., Levshakov S. A., et al. 2006, A&A, 
454, 151 

Wolfe A. M., Gawiser E., Prochaska J. X. 2005, ARA&A, 43, 
861 



© 2012 RAS, MNRAS 000,[THTT] 



